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The shape and the inverse participation ratio (IPR) of lo- 
cal spectral density (LSD) are studied for a generic isolated 
system of coupled quantum states, the Hamiltonian of which 
is represented by a band random matrix with the disordered 
leading diagonal. We find for the matrices with arbitrary 
small band that the lack of ergodicity for LSD can be associ- 
ated with an exponential increase in IPR with the ratio v/A c 
(v - the root of mean square for off-diagonal matrix elements, 
A c - the energy spacing between directly coupled basis states). 
Criterions specifying transition to localization and ergodicity 
for LSD are considered. 
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I. INTRODUCTION 



The statistical properties of quantum system attract 
considerable attention in the broad field of modern 
physics ranging from nuclear, atomic and molecular 
physics to condensed matter and quantum computing. 
Of a special interest is the properties specifying the quan- 
tum chaos border for a system of many-body quantum 
states. Extensive investigations of many-body interact- 
ing systems such as nuclei, many-electron atoms, quan- 
tum dots, quantum spin glasses and quantum computer 
models have shown that the border is associated 

with a crossover of level spacing statistics from the Pois- 
son distribution to Wigner-Dyson one and with a transi- 
tion to ergodic distribution for the eigenstate wavefunc- 
tion over the large number of basis levels as in Random 
Matrix Theory (RMT) (lO0. In a sense the interac- 
tion leads to dynamic thermalization without coupling 
to an external thermal bath. According to [p|,|6||T^l this 
crossover takes place when the coupling matrix elements 
of the Hamiltonian become comparable to the energy 
spacing between directly coupled states. 

We address this problem to a generic system of coupled 
quantum states, the Hamiltonian matrix H of that in- 
cludes the leading diagonal with disordered random val- 
ues and the random off-diagonal elements inside the band 
of size b. This band random matrix with the disordered 
diagonal (BRMDD) is a reasonable model describing the 
systems with strong imperfection in the basis state ener- 
gies. Models based on BRMDD were applied to study the 
electron transport problem p^Jl4] | and the problem of in- 
teracting particles in a random potential ]H)|-|l7|] . The re- 
sults obtained for these models arc of an obvious interest 



also in analyzing such few-freedoms physical objects as 
the vibrational quasicontinuum of polyatomic molecule 
p8[ and a quantum computer 

The statistical properties of BRMDD-based systems 
were studied in some details |I1|pi||l9| p4| . The shape, 
the localization length and the inverse participation ratio 
(IPR) for eigenfunction have been investigated with the 
help of numerical simulations 15 1(||22| and the super- 
symmetry approach ^,^3|. These investigations have 
exhibited the Lorentzian shape for local spectral den- 
sity (LSD) in circumstances when a non-perturbativc 
localization regime is realized. The state-state interac- 
tion strength, at which the eigenstates are extended over 
the whole matrix size N and the eigenenergy level spac- 
ing statistics has the Wigner-Dyson form, have been re- 
vealed for BRMDD with sufficiently large band (when 
26 + 1 > VN) @^||). Unfortunately the ergodic 
properties of LSD have not been investigated in details. 
An open question is also the quantum chaos border for 
BRMDD with arbitrary small band. 

In the paper we study localization and ergodicity prop- 
erties of LSD to be obtained from BRMDD with arbitrary 
small band. The LSD was introduced in 1955 by Wigner 
p5j and successfully employed in RMT to describe sta- 
tistically the localization effects for complex quantum 
systems |H],|ll| (including the systems, represented by 
band random matrices with the reordered leading diag- 
onal [^6 27 1). This quantity is the Fourier transform for 
the correlation function C(t) = (0| exp(— iHt/h) |0) and 
specifies spreading of the energy, initially concentrated in 
a specific unperturbed basis state |0), between the eigen- 
states due to state-state interaction. Generally the LSD 
may be characterized by the width T measuring the en- 
ergy scale, in that the individual state |0) is localized. 
The number of eigenstates populating this scale is given 
in terms of the product 



£,e = PET 



(1) 



designated here as an ergodic localization length of LSD 
(Pe is the eigenstate density). Hence the quantity £ e 
specifies the greatest possible number of eigenstates, 
where the basis state |0) can be effectively admixed. The 
ergodic properties for the system can be identified with 
the structure of LSD. The non-ergodic LSD is a strongly 
fluctuating spiked function and the IPR £ipr, which 
gives the actual number of eigenstates involving the state 
|0), is low in comparison with £ e - In the ergodicity case 
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the LSD is monotonic and the number £,ipr approaches 
the value of £ e . 

The main attention of our study is paid to the per- 
turbative and localized regimes for LSD. Therefore we 
restrict the study to a moderate strength of state-state 
interaction. This restriction implies that the width T is 
essentially small as compared to the energy scale bound- 
ing the location of the eigenstate levels (£,ipR,£, e « N). 
The finite-size effects are ignored. Hence the quantities 
£ e and £jp_r can be considered as functions of the inter- 
action strength and the relative width of band. 

We investigate the LSD shape and determine the er- 
godic localization length. Then we analyze the behav- 
ior of IPR for ergodic and non-ergodic LSD. From this 
analysis we obtain criterions specifying the transition to 
localized and ergodic regimes for LSD. 



When the coupling matrix elements are weak and 
v/A c <C 1 the perturbative regime is realized. At the 
regime the Hamiltonian (||) can be treated within the 
framework of perturbation theory. The condition 1 <C 
v/A c <C Vb features the regime, at that the eigenstatcs 
are localized. The strong interaction at v/A c \/b is 
responsible for the delocalizcd regime for cigenfunction. 
In the regime the eigenstates are extended over the whole 
matrix size N. 

The eigenfunctions \a n ) and the eigenenergies E n are 
obtained from diagonalization of (J2|) . The quantity 
W n o = | (a n 1 0)| 2 gives a probability to find the probing 
state |0) in the eigenstate \a n ). The numerical simulation 
is performed for the wide range of matrix parameters: 
2 < N < 3200, 1 < b < 1600 and 10~ 3 < v/A < 10 2 . 
The number of disorder realizations lies in the range from 
100 to 1000. 



II. MODEL DESCRIPTION 

We consider a BRMDD, that represents the Hamil- 
tonian matrix H in the basis of unperturbed states \k) 
(k = -K, . . . , AT) for isolated system of N = 2K + 1 
quantum states. The BRMDD is a real symmetric matrix 
with statistically independent random elements 



H„ 



5-mk + V m 



(2) 



where off-diagonal elements V m k = Vkm specify the state- 
state interaction. The values of V m k are distributed 
uniformly in the interval [— V, V] with (V m k) = and 



„2 _ 



(V 2 



V /3 if \m — k\ < b or are zero otherwise. 



The diagonal elements E^ corresponding to energy lev- 
els for the states \k) are uniformly distributed according 
to the Poisson statistics with the mean spacing A be- 
tween adjacent energy levels: —KA < E?' < KA. The 
energy level E^ for a probing state |0) is located in the 



(0)\ 



E, 



(0) 



0. 



midpoint of interval [—KA,KA\: (E i 

The perturbation of an individual basis state depends 
on the relative level position for unperturbed states, 
those are coupled directly with the state by interaction. 
Therefore the strength of state-state interaction for BR- 
MDD can be specified in terms of the ratio v/A c , where 



A 

1 



(3) 



is the energy spacing between directly coupled unper- 
turbed states. Here the relative band width 



b 

K 



(4) 



gives the relative number of states coupled directly with 
the individual state. 

Depending on the relation between v/A c and b, three 
important regimes can be distinguished plJT^,p2l,p3| . 



III. SHAPE OF LSD 

The shape of LSD is studied for the eigenenergy scale, 
which is no larger than 5r. For this scale the interaction- 
induced variations in level density are negligible and 
the eigenenergies are homogeneously distributed with the 
density pe ~ 1/A. Then the LSD can be defined as 



Pw {E) = pe 



(EnWnoSjE - E n )) 

(E n 6(E-E n )) ■ 



(5) 



where (...) means the averaging over disorder (that is, 
over many random matrices). The numerical investi- 
gation confirms that the LSD is described by the well- 
known Breit-Wigner distribution of the Lorentzian shape 



Pbw(E) 



F 



1 



2tt E 2 +T 2 /4 



(G) 



(see the inset in Fig.|l|) with the width T. This shape is 
observed for any v, A c , b and satisfying £ e <§; N (in- 
cluding the parameter region where r <C A). We find the 
ergodic localization length £ e from (|l]) where the width 
r is estimated by fitting of the Lorentzian (Q) to an av- 
eraged LSD. Our analysis of numerical data obtained for 
different values of v, A c and (3 shows that the quantity £ e 
can be considered as a function of the single parameter 



(7) 



A c vr 



which defines essentially the localization of LSD in the 
eigenenergy scale. When the coupling of the basis states 
by perturbation is sufficiently strong and q » 1 the 
quantity £ e satisfies the Fermi golden rule: £ e A = T m 
2-rtq 2 A = 2ttv 2 / A c > A. At low values of q (q « 1) 
the length £ e can be described by a linear function of q: 
£ e w 2q. The width T is then much smaller as compared 
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to the spacing A. Fig.[j] demonstrates the behavior of £ e 
in the wide range of q at three fixed magnitudes of (3. 
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FIG. 1. Dependence of the ergodic localization length £ e 
on q for j3 =1 (+), 0.5 (O) and 0.04 (A). The solid line gives 
the approximation (9). The inset shows a semilog plot for 
the LSD p w {E) for (3 = 1 at v/A c = 0.089 (+, 1000 disorder 
realizations), v/A c = 0.45 (O>100 disorder realizations) and 
v/A c — 2.82 (A, 100 disorder realizations). The solid lines 
give the Breit-Wigner distribution (6) with £ e =0.05, 1.3 and 
50. 

The linear dependence of £ e and T on the strength of 
state-state coupling for small q can be explained in terms 
of a simple perturbation approach to the cigcncncrgy lcv- 



E n — E, 



(0) 



< 



V < 



els E n located in the vicinity of Eq 

A. In the approach the eigenstate \a n ) is considered as 
a superposition of the nearest basis states |0) and |n). 
The probability W n o to find the basis state |0) in the 
eigenstate \a n ) can be then estimated from a two-state 
(the states |0) and \n)) approximation. According to this 
approximation the dependence of W n o on E n for an in- 
dividual disorder realization is Lorcntzian-likc 



W n0 



V 2 



V* 



E n — E. 



(0) 



(8) 



with the width of 2 |V^q|. If the states |0) and \n) are 
not coupled by interaction the dependence (||) is con- 
sidered as a function, the width of that is negligible in 
comparison to v. The probability that the states |0) and 
\n) are coupled by interaction is equal to (3. As a result 
the shape of disorder-averaged LSD is described by the 
Lorentzian-like contour with the width T f=a 2j3v. It seems 
that this perturbation approach can be applied also for 



the eigenstates with E n — E { 



P (0) 



A, when the proba- 



bility Wno is essentially small as compared to the contri- 
butions of others basis states. In this limit Eq. (^) gives 
Wno ~ V 2 / (En —Eq) 2 , as predicted in the framework 
of a standard perturbation approach for non-degenerate 
states (see, for instance, in ps{]). 

In order to estimate the localization properties of LSD 
at intermediate values of q we approximate the behavior 
of ergodic length for the wide range of q (10~ 3 < q < 15). 
The analysis of numerical data shows that the length £ e 
can be well approximated by 



& « L iqy /1 + (L 2 q) 2 



(9) 



with the fitted coefficients L\ = 2.01 w 2 and L 2 — 
3.16 » 7T. According to (||) an accurate boundary be- 
tween the linear and quadratic dependencies of £ e on q 
is obtained at q w I/L2 where £ e w 0.9. Notice that this 
boundary separates the ranges where the value of £ e is 
low or high as compared to 1. 



IV. PARTICIPATION RATIO FOR LSD 

The IPR i IPR = ((£»! W§ I 2 ),)" 1 is associated 
with the parameters q and (3, which specify correspond- 
ingly the ergodic localization length £ e for LSD and the 
band width for BRMDD. This association depends on 
whether or not the LSD is ergodic. In the ergodic case the 
probing state is monotonically spread over eigenstates in 
the energy scale T. Then the IPR is a function of the only 
parameter q and approaches the ergodic length £ e . The 
ergodic LSD can be realized from BRMDD with large 
band. Our analysis of calculated data for the matrices 
with (3 ~ 1 shows that the best approximation for ergodic 
IPR is 



6p fl w 1 + Dtqy/l + (D 2 q) 2 a 1 + £ e 



(10) 



The fitted coefficients Di = 3.16 and D 2 = 1.94 are 
obtained for full Hamiltonian matrices (0) with = 1 
where all 2K basis states are directly linked with the 
state |0). Notice that the approximation (10) gives the 



maximum allowable magnitude of IPR for BRMDD at 
fixed q. 

An essentially non-ergodic LSD can be obtained from 
BRMDD with /3 < 1. Then the IPR is significantly low 
in comparison to £ e and should be considered as a func- 
tion depending both on the parameter q and on (3. The 
analysis of numerical data for matrices with low magni- 
tudes of q and (3 shows that the lack of ergodicity for LSD 
can be clearly identified with an exponential increase in 
IPR with the quantity q\f]3 = v/A c 



(11) 



£ipr ~ exp Cgv/3 



At fixed j3 the IPR is shown in Fig.^ to increase expo- 
nentially with the rate C\f]3 (C a 3.0 -j- 3.2) as the pa- 
rameter q rises if £,ipr <C £ e . The increase is observed up 
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to the point where £,ipr ~ £ e ■ Then the IPR approaches 
asymptotically the length £ e from below with increasing 
q. At high magnitudes of q and j3, when the exponential 
(O) gives higher values of £,ipr in comparison to £ e , the 



LSD is ergodic and the IPR can be determined from (10) 




FIG. 2. Dependence of the IPR £ IPR on q for /3 =1 (full 
circle), 1/2 (O), 1/4 (+), 1/8 (A), 1/16 (full triangle) and 
1/32 (□). The dashed lines give the exponential (11) for 
/3 = 1/4 (C = 2.99), 1/8 (C = 3.15), 1/16 (C = 3.19) and 
1/32 (C = 3.13). The solid line shows the ergodic approxi- 
mation (6) of IPR. 




FIG. 3. Dependence of the IPR £,ipr on the product qyf]3 
for the broad range of q (0.01 < q < 10) and /?(/?> 0.001) sat- 
isfying the requirement £ipr < £ e /2.7 (+). Solid line shows 
the approximating exponential (11) with C = 3.15. 



For the broad range of q and (3 satisfying to £,ipr < 
£ e /2.7 the mean square fitting gives the coefficient C — 
3.15±.01. FigJ| demonstrates a good accordance between 
the calculated and fitted values of IPR for non-ergodic 
LSD. Unfortunately we can give no heuristic explanation 
of the exponential rise in IPR with q\f]3. Notice only that 
this rise is observed under conditions when the eigen- 



states are essentially localized and l s b/N ft 
(l s b is the localization length for eigenstate). 



q 2 p 2 « 1 



V. CRITERIONS OF LOCALIZATION AND 
ERGODICITY OF LSD 

We associate the localized regime for LSD with a high 
magnitude of IPR. For BRMDD with arbitrary small 
band the requirement £ipr 3> 1 is seen from ( |Io|) and 
( |i~T| ) to result in 



v/A c = q^fp > q d y% = 1/Co ~ 1/tt, 



(12) 



no matter what the ergodic properties of LSD and the 
magnitude of (3. The parameters qd and (3d defines then a 
localization border on the plane {q, 0), that is a boundary 
between regions where the perturbative {£,ipr ~ 1) or 
localized (£ipr 3> 1) regimes are realized. Notice that 
as for LSD, the criterion ( [l2| ) specifies a transition to the 
localized regime for eigenstate 15 l?j, ^2] , p3| . 

The attainment of ergodic LSD for BRMDD is defined 
by the relative width of band. For the matrices with (3 ~ 
1 the localized LSD is nearly ergodic and the ergodicity 
criterion can be represented by the condition (|l2]). This 
condition is in good agreement with a quantum chaos 
border for many-body systems and gives the crossover 
of level spacing statistics from Poisson to Wigner-Dyson 
distribution for a nucleus model (Co = 3.0 and for 
a generic model of quantum computer (Co — 2.5 0). In 
( |j~2| ) the factor Co is assumed to be equal to C « D\ w it. 

For BRMDD with j3 <C 1 we should take into account 
the exponential growth ( ]i~l| ) of IPR with v/A c . The cri- 
terion of LSD ergodicity can be determined then from 
the requirement exp(Cgv / /3) > £e, that results in: 



v/A c = q^fp > q eA //% = 



In (2Trq 2 ) _ In (2nql 



C 



(13) 



where the values of q e and f3 e specify an ergodicity border 
on the plane (q,f3). Notice that the condition ( |i~3| ) is a 
more severe restriction imposed on the parameters q and 
(3 than the criterion (p^|). 

Our study reveals some distinctive features in behavior 
of LSD for BRMDD with arbitrary small band. These 
features are outlined on the plane (q, f3) in Fig.|| to de- 
pend on the relation between v/ A c and (3. As for eigen- 
state, the perturbative regime for LSD can be associated 
with the requirement v/A c <g; 1/3. In contrast to eigen- 
state, the localization properties of LSD is essentially de- 
fined by the magnitude of q. In the limit q <C 1/3 (that 
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is, v/A c <§C \f]3/2>) the ergodic localization length is a lin- 
ear function of q and £ e w 2v/A c ,/j3. Due to weak state- 
state interaction the LSD is a delta-like function confined 
in an eigenstate: £ e -C £ipr ~ 1. At higher values of y 
(1 < 3y < or y7? < 3v/A c < 1) the quantity £ e 

satisfies the Fermi golden rule: £ e » 2ttv 2 / A 2 f3 . Despite 
the wide spreading of basis states over the eigenenergy 
scale, the actual number of eigenstates associated with an 
individual basis state for this parameter region remains 
small: £ e > £,ipr ~ 1. 



10" 



tribution (see Fig.||) to be associated usually with quan- 
tum chaos. In a sense the condition ( |l3| ) can be consid- 
ered as a quantum chaos criterion. 
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FIG. 4. Localization regimes depending on the parameters 
q and /3. The fat dashed line shows the localization border 
(12). The fat solid curve gives the LSD ergodicity border (13). 
The dotted line indicates a boundary between the regions with 
£ e < 1 and £ e > 1. The points A (g = 5.55, = 0.016) and B 
(q — 5.55, /3 = 0.125) give the parameter sets used to obtain 
the level spacing statistics represented in Fig. 5. 

The criterion (|l2| ) specifies a transition from the per- 
turbative to localized regime both for eigenstate and for 
LSD. Depending on the ergodicity properties of local- 
ized LSD, one can distinguish two important parameter 
ranges. An area, located between the borders (|lj) and 
( |l3| ) on the plane (q,f3), corresponds to the parameter 
region 1 < 3v/A c < ln(£ e ) where £ e « 2nv 2 /A 2 c f3 » 
£iPR « exp(Cu/A c ) >> 1. For this region the neighbor- 
ing eigenstates are weakly coupled by interaction since 
the density £ipp/T of eigenstates, those involve the prob- 
ing state, is negligible in comparison to the total state 
density pe — £e/r. It means that the level spacing 
statistics P{S) should exhibit the Poisson-like distribu- 
tion as it shown in Figj^. We identify this region with 
a non-ergodic localized regime for LSD. The regime of 
ergodic localization for LSD at v/A c ^> ln(£ e )/3 is rep- 
resented by a domain above the ergodicity border (|l3| ) 
where £ e « 2ttv 2 / A 2 (3 ~ £,ipr >> 1- For this domain the 
level spacing statistics is close to the Wigner-Dyson dis- 
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FIG. 5. Level spacing statistics P(S) at £ e = 193.5 for 
non-ergodic (/3 = 0.016, £ipr = 10.3, solid line stairs) and 
ergodic (f3 — 0.125, £ipr = 150.5, fat line stairs) LSD lo- 
calization regimes, represented correspondingly by the points 
A and B in Fig. 4. The dashed curves show the Poisson and 
Wigner-Dyson distributions. 



VI. CONCLUSION 

We have analyzed the properties of LSD for a generic 
conservative system of coupled states, the Hamiltonian 
of which is represented in terms of a BRMDD. As for 
eigenstate, the transition from the perturbative to local- 
ized regime for LSD can be associated with the same 
requirement to be imposed on the strength of state-state 
interaction. Nevertheless, our study has revealed some 
distinctive features in the properties of LSD. We have 
determined localization and ergodicity borders for LSD 
obtained from BRMDD with arbitrary small band. For 
the matrices with large band (/? ~ 1) the borders co- 
incide with each other and are in good agreement with 
the quantum chaos border for some generic many-body 
systems. For the matrices with small band (f3 <C 1) the 
criterion of LSD ergodicity is a more severe requirement 
imposed on the band size and the interaction strength 
than the localization criterion. In the parameter range 
between the borders the IPR is found to increase expo- 
nentially with the ratio v / A c , which specifies the strength 
of state-state interaction. We suppose that such an effect 
can be observed for any isolated system of coupled states 
with strong imperfection in non-perturbed energies. 

Author thanks D. Shepelyansky and V. Churakov for 
useful discussions and remarks. 
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